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Abstract 

P odd anapole moment of the deuteron is found in the chiral limit, 
m,r — ^ 0. The contact current generated by the P odd pion exchange 
does not contribute to the deuteron anapole. Being combined with 
usual radiative corrections to the weak electron - deuteron interaction, 
our calculation results in a sufficiently accurate theoretical prediction 
for the corresponding effective constant C2d- The experimental mea- 
surement of this constant would give valuable information on the P 
odd ttNN constant and on the s-quark content of nucleons. We cal- 
culate also in the same limit m^r — > the deuteron P odd and T odd 
multipoles: electric dipole moment and magnetic quadrupole moment. 
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1 Introduction 

The investigations of nuclear parity violation, both theoretical and experi- 
mental, have already a long history. New light on this problem is shed by 
observation of the nuclear anapole moment (AM) of ^^^Cs in atomic experi- 
ment [0] . The result of this experiment is in a reasonable quantitative agree- 
ment with the theoretical predictions, starting with [Q, ^, if the so-called 
"best values" |Q are chosen for the parameters of P odd nuclear forces. 

The AM is a rather peculiar multipole in the following sense (for a more 
detailed discussion see, for instance, 0). The interaction of a charged probe 
particle with an anapole moment is of a contact nature. Therefore, for in- 
stance, the interaction of the electron with the nucleon AM, being on the 
order of aG, cannot be distinguished in general case from other electromag- 
netic radiative corrections to the weak electron-nucleon interaction due to the 
neutral currents. And in a gauge theory of electroweak interactions only the 
total scattering amplitude, i.e., the sum of all diagrams on the order of aG, 
is gauge- invariant, independent of the gauge choice for the Green's functions 
of heavy vector bosons (here a is the fine-structure constant, G is the Fermi 
weak interaction constant). No wonder that, generally speaking, the AM 
of an elementary particle or a nucleus is not gauge- invariant, i.e., physically 
well-defined, quantity. However, there is a special case where the anapole 
moment has a real independent physical meaning. In heavy nuclei, of course 
^^^Cs included, the AM is enhanced ~ A^/^ [A is the atomic number), as 
distinct from common radiative corrections. By the way, it means that there 
is an intrinsic limit for the relative accuracy, ~ A"^/^, with which the AM 
of a heavy nucleus can be defined at all. For ^^^Cs this limiting accuracy is 
about 4%. 

There is one more object, the deuteron, whose anapole moment could 
make sense for a sufficiently large P odd vrNN constant ||^. The deuteron 
is a loosely bound system of a relatively simple structure. Therefore, there 
are all the reasons to believe that its AM is induced mainly by the P odd 
TT-meson exchange, pion being the lightest possible mediator of the nucleon- 
nucleon weak interaction. The problem of the deuteron AM was discussed 
phenomenologically in [2,6-8]. In the present work the deuteron AM, as 
induced by the P odd vr-meson exchange, is explicitly expressed through the 
P odd vrNN coupling constant. The same problem was considered in recent 
paper 0. The result of the original version of was much smaller than 
ours because a leading contribution, that of the isovector magnetic moment 
of the nucleon, was omitted in it. After acquaintance with the preprint of 
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the present work, the authors of [^] corrected their result (see their revised 
preprint 0, Section VI. Erratum and Addendum). On the other hand, under 
the influence of 0, we have refined our own calculations with the results 
presented below. 

The obtained result for the deuteron AM is singular as l/m.,^ in the limit 
m^r — > (of course, when going over to this limit, one should keep the 
deuteron radius larger than the Compton wave length of the pion). Being 
combined with the radiative corrections to the weak electron - deuteron 



scattering amplitude [jTO[, which are regular in m^, our calculations result in 
a sufficiently accurate value for the corresponding effective constant C2d- 

We also calculate here the P odd and T odd electromagnetic moments of 
the deuteron. 



2 The deuteron anapole moment 

It is convenient to start the discussion with the nucleon AM in the chiral limit. 
It was shown in 1980 by A.I. Vainshtein and one of the authors (I.Kh.) to be 
given in this limit by the diagrams 1 and 2. The circle on the nucleon lines 




— e X — — X e — 

p(n) n(p) p(n) p (n) n (p) p (n) 

Fig.l Fig. 2 

refers to the usual strong interaction vrNN vertex (coupling constant g\/2), 
the cross describes the P odd weak vrNN interaction (coupling constant g\/2). 
The result for the nucleon AM is 

egg ( Q mp\ 
a^r = ap = a„ = - — 1 In cr. (1) 

llniprnT^ \ tt nip m.,^ J 

The diagrams discussed lead to the same result for a proton and neutron 
since under the permutation p ^ n the strong coupling constant g does not 
change, and the weak one g changes sign together with the charge e of the tt- 
meson (we assume e > 0, exact definitions of the strong and weak interaction 
Lagrangians and coupling constants g and g are given below). Being the only 
contribution to the nucleon AM, which is singular in m^r, the result (|I|) is 
gauge-invariant. In this respect, it has a physical meaning. 
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Unfortunately, in spite of the singularity in m^, the corresponding con- 
tribution to the electron-nucleon scattering amplitude is small numerically 
as compared to other radiative corrections to the weak scattering amplitude. 
Indeed, the radiative corrections to the effective constants C2p^n of the proton 
and neutron axial neutral-current operators G/v^ C2p,nCrp,n are [^] 

C^p = 0.032 ± 0.030 , = -0.018 ±0.030. (2) 

In the same units G/\/2, the effective axial constants induced by the elec- 
tromagnetic interaction with the proton and neutron anapole moments (|ip, 
is 

C^^^^ = -aaN{\e\G/V2)-'^ = 0.07 x 10^^. 

At the "best value" g = 3.3 x 10^^ (strongly supported by the experimental 
result for the ^^^Cs anapole moment) we obtain 

C;,n = 0.002. (3) 

With this value being much less than both central points and error bars in 
(0), the notion of the nucleon AM practically has no physical meaning. This 
is why the result ([l|) was never published by the authors. It is quoted in 
book [^] (without the logarithmic term) just as a theoretical curiosity. The 



logarithmic term in the nucleon AM is discussed in ||TT . 

However, the situation with the deuteron AM is quite different. Not only 
the proton and neutron AMs add up here. The isovector part of the radiative 
corrections is much smaller than the individual contributions and Cj^, 



and is calculated with much better accuracy ||1(| : 

Cm = C2p + = 0.014 ± 0.003 . (4) 

Moreover, there is the already mentioned, qualitatively new contribution, 
due to the isovector magnetic moment of the nucleon, which dominates nu- 
merically the deuteron AM. Thus acquires a real physical meaning. 

Let us go over now to the problem itself. The Lagrangians of the strong 
ttNN interaction and of the weak P odd one, Lg and L^, respectively, are 
well-known: 

Ls = g[\/2 {pi-f^n 7r+ ± nij^p vr") ± ( pijr,p - nij^^n) 7r°] ; (5) 
Ly^ = g V2i {pnn^ — npn^). (6) 



Our convention for 75 is 




( , ; (7) 



the relation between our P odd ttNN constant g and the common one /i^^rAf 
is ^^2 = /iS^. 
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The effective nonrelativistic Hamiltonian of the P odd nucleon-nucleon 
interaction due to the pion exchange is in the momentum representation 



F(q) = ^ J^^i^Nln.N,) {N^T.+N,). (8) 
nip ml + 



Here 



1 

2' 

is the deuteron spin, q = - = -(p^ - pg) = p^ - p^ = -(p^ - p„). 
Let us note that the P odd interaction (H) which interchanges the proton and 
neutron, when apphed to the initial state aj,(ri)ajj(r2)|0) transforms it into 
ajj(ri)aj,(r2)|0) = — aj^(r2)ajj(ri)|0). On the other hand, the coordinate wave 
function of the admixed ^Pi state is proportional to the relative coordinate 
r, which we define as — r„. Therefore, it also changes sign under the 
permutation p ^ n. Thus, for the deuteron the P odd potential can be 
written in the coordinate representation as a simple function of r = rp — r„ 
without any indication of the isotopic variables: 

V{r) = . V) '.^^tUhlll (9) 

The above expressions are rather standard. As standard is our sign conven- 
tion for the coupling constants: g = 13.45, and ^ > for the range of values 
discussed in [If. 

The discussed P odd interaction V generates a contact current j*^. To 
obtain the explicit expression for it, we have to consider V in the presence 
of the electromagnetic field. Its including modifies the proton momentum: 
p ^ p — eA, which results in the shift q ^ q + eA in the interaction (H). 
Then in the momentum representation the contact current is 

. , 9V(q) d 2gg I(q + eA) 

J (q 



OA OA nip nil + (q + ^A) 
2e^ \_}: 2q(Iq) ] 



TUp y ml + q^ [ml + q^ 

In the last expression we have neglected the dependence of the contact current 
on A. In the coordinate representation it equals 

m = ^riivf-^. (11) 
lump r 

Let us derive at first a general structure of the deuteron AM generated by 
a P odd np interaction, assuming only that the deuteron is a pure ^5*1 state, 
bound by a spherically symmetric potential. We follow here essentially the 
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line of reasoning applied in (see also book to the problem of a single 
proton in a spherically symmetric potential. In this case the formula for the 
AM operator is [Q] 

a = — {/ipr X o- - ^ [1^, r]} + ^ r X [r X j'^], (12) 
nip 3 3 

with the proton magnetic moment fip = 2.79. In the case of the deuteron 
this formula generalizes to 

ad = ^ {r X (/ipo-p - ixno-n) - - [l^ r]} + - r X [r X j^], (13) 
fin = —1.91 is the neutron magnetic moment. Both AM operators ([T2| ) 



and (13) are orthogonal to r (neither of them commutes with r, so the or- 
thogonality means here that ar + ra = 0). Therefore, the contact current 
(p!T|) generated by the P odd pion exchange and directed along r, does not 
contribute to the nuclear AM. 

Let us present now the wave function of the deuteron ^5*1 state as V'o(^)X) 
where x is the spin wave function for 1 = 1 (we neglect here and below a 
small ^Di admixture in the deuteron). If the P odd interaction conserves 
the total spin I of the deuteron, the ^Pi state admixed by it can be written 
as i{lr/r)ilji{r) (both radial wave functions, ipo{r) and ipi{r), are spherically 
symmetric). Simple calculations demonstrate that the deuteron AM, as in- 
duced by the operator (|T3|), is in the absence of the contact contribution 



ad = ^ (^/ip ~ ~ ^) / drriJo{r)iJi{r). (14) 



So, under the assumptions made, the deuteron AM should depend on the 
universal combination (/ip — /i„ — 1/3). 

We confine mainly in our calculation to the naive zero-range approxima- 
tion (ZRA) for the deuteron wave function: 

/(o)/ N exp (-/tr) 

^o'ir) = ^— . (15) 

Here n = ^vHpS] e = 2.23 MeV is the deuteron binding energy. 

The P odd correction to the deuteron wave function due to V{r) will be 
found in the common stationary perturbation theory. In the same ZRA the 
admixed ^Pi states of the continuous spectrum are free. Moreover, we can 
choose plane waves as the intermediate states since the perturbation V{r) 
selects by itself the P-state from the plane wave. Thus obtained first-order 
correction to the wave function is 

*'C^) - / (5^ rF^k /*'»-*''''(--')*»(^'). (16) 
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Rather lengthy calculation leads to the following expression for the matrix 
element of the radius- vector: 

driljr)(r)rilji(r) = - — ^ttt, (17) 

where ^ = n/mT^ = 0.32. With this matrix element and the operator ([T3|) 
one obtains easily the following result for the deuteron AM: 

in accordance with the general formula (|14|)- Our overall factor at the 
structure {fip — /i„ — 1/3) is the same as that at (/Xp — fin) in the re- 
vised version of 0. However, the corresponding total result obtained in 0, 
even in it revised version, is not proportional to the universal combination 

(yUp - jJn- 1/3). 

In fact, the range 1/m^ of the P odd interaction (||) is quite comparable 
to the range of the usual nuclear forces. Therefore, it is, strictly speak- 
ing, inconsistent to use the zero-range approximation for calculating effects 
induced by the perturbation V^(r). Still, numerical estimates made with a 
model deuteron wave function which has somewhat more realistic properties, 
indicate that the error introduced by using the ZRA does not exceed 20%. 
As to other sources of P violation, different from the pion exchange, there 
are no reasons to expect that in the case of deuteron their neglect creates a 
serious error if the P odd tcNN coupling constant g is at least comparable 
to its "best value". 

It looks reasonable to combine the potential contribution (|TE|) with the 
additive contribution of the nucleon anapole moments, which according to 

aJ' = a, + a„ = --i»^fl-«^!^;Llnn^Ell. (19) 



In this way we arrive at the final result for the deuteron AM in the chiral 
limit: 



egg 

ad = - 



0.49 (/xp - /in - ^) + 0.46 



I = -2.60 I. (20) 



This result includes all contributions to the P odd amplitude of ed-scattering, 
which are singular in m^r, and thus is gauge-invariant, independent of the 
gauge choice for the Green's functions of heavy vector bosons. 

Finally, let us compare the contribution of ( pUD to the P odd ed scattering 
amplitude which is due to the usual radiative corrections, nonsingular in m„. 
For the deuteron the axial operator looks as follows: 
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The contributions to the isoscalar axial constant C2d originate from the 
anapole moment, from usual radiative corrections nonsingular in m.,^, and 
from the admixture of strange quarks in nucleons . The magnitude of the 
s-quarks contribution is extremely interesting, but highly uncertain. As to 
the usual radiative corrections, their contribution to this constant is found 
in [|10| with good accuracy (see (^): = 0.014 ± 0.003. In the same units 
G/ a/2 the effective axial constant induced by the electromagnetic interaction 
with the deuteron AM (^) is 

C^^ = aad{eG/V2)-^ = 0.44 x 10^^. (21) 

At the "best value" g = 3.3 x 10^^ (strongly supported by the experimental 
result for the ^'^^Cs anapole moment) we obtain 

C^a = 0.014 ± 0.003. (22) 

We use here the above estimate of 20% for the accuracy of our calculation 
(for given g). The numbers in (^) and (^2]) are quite comparable, and taken 
together result in the following value of the total effective constant: 

C^d = Cl^ + C^^ = 0.028 ± 0.005. (23) 

We are fully aware of the extreme difficulty of the experimental mea- 
surement of the constant C2d- However, with such a good accuracy of the 
theoretical prediction (p3|) , this experiment becomes a source of the valuable 
information on the P odd vrNN constant and on the s-quark content of nucle- 
ons. As it was 15 years ago, now again "C2d seems to be the most interesting 



parity-violating parameter accessible to atomic-physics experiments" |T0 
although by rather different reasons. 

Of course, if necessary the accuracy of our prediction ( p2D can be im- 
proved by using a more detailed and realistic description of the deuteron. 
On the other hand, the accuracy of radiative corrections (^ can be also im- 
proved, at least by using much more precise modern experimental values of 
the parameters of the electroweak theory. 



3 The deuteron P odd, T odd moments 

The problem of the deuteron P odd, T odd multipoles: electric dipole, mag- 
netic quadrupole, and the so-called Schiff moment, was treated phenomeno- 
logically in |l^. Now we will calculate the electric dipole and magnetic 



quadrupole moments within the approach applied above to the anapole. As 
to the Schiff moment, strong cancellations occur when calculating its value 



for the deuteron [14]. Therefore, one cannot expect reasonable accuracy for 
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it with our ZRA deuteron wave function, and we will not consider here this 
problem. 

As distinct from the P odd, T even interaction, there are three inde- 
pendent P odd, T odd effective vrNN Lagrangians. They are conveniently 
classified by their isotopic properties: 

(24) 
(25) 



AT 


= 0. 




= 90 


\pl {pn TT^ 


+ npn ) + {pp — nn) ti^] ; 


|AT| 


= 1. 




= 91 


NNtt^ = 


gi {pp + nn)TT% 


|AT| 


= 2. 




= 92 


{NtNtz 


- 3Nt^Ntt^) 








= 92 


[ \pl {pn TT 


^ + npn^) — 2{pp — nn ) 



(26) 

Since the possible values of the isotopic spin for two nucleons is T = 0, 1 
only, the last interaction, with |AT| = 2, is not operative in our approach. 

The effective P odd, T odd proton - neutron interaction is derived in the 
same way as in the AM problem. In the momentum representation it looks 
as follows: 

= ^ 2I , [{^9o-9i)<T, - iSgo + 9i)<Tn]. (27) 
In the coordinate representation it is 

W{r) = -^[{3go-9i)^p - ( 3^o + ) J V . (28) 

Simip r 

The calculation of the deuteron EDM d^, i.e., of the er^ = er/2 expec- 
tation value, goes along the same lines as that for the anapole moment and 
results in 

d=-^^-l±i-I. (29) 
127rm^ (1 + 20^ ^ ^ 

The magnetic quadrupole moment (MQM) operator is expressed through 
the current density j as follows (see, for instance, 0, p!5|] ): 

This expression transforms to 



e 
2m 



js/i r^cTn + r^cxm - ^ (err) + 2g(r„/„ + r„/m)| . (31) 



Here /x is the total magnetic moment of the particle, q is its charge in the 
units of e. The magnetic quadrupole moment is the expectation value M. of 
the operator M^z in the state with the maximum total angular momentum 
projection = /. 

In our case, due to the spherical symmetry of the deuteron nonperturbed 
wave function, the orbital contribution to Mmn vanishes. The contact current 
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generated by the P and T odd charged pion exchange, here is also directed 
along r, and thus does not contribute to MQM. So, the deuteron magnetic 
quadrupole moment originates from the spin term in (^). It equals 

M = - ./XI 2 [ ( 3^0 + ^1 ) + ( 3^0 - ^1 ) ) ]. (32) 
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